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We show that ultracold atoms can be controlled in multi- 
band optical lattices through spatially periodic Raman pulses 
for investigation of a class of strongly correlated physics re- 
lated to the Kondo problem. The underlying dynamics of 
this system is described by a spin-dependent fermionic or 
bosonic Kondo-Hubbard lattice model even if we have only 
spin-independent atomic collision interaction. We solve the 
bosonic Kondo-Hubbard lattice model through a mean-field 
approximation, and the result shows a clear phase transi- 
tion from the ferromagnetic superfluid to the Kondo-signet 
insulator at the integer filling. 

PACS numbers: 03.75.Fi, 03.67.-a, 42.50.-p, 73.43.-f 
Ultracold atoms in optical lattices have recently re- 
ceived a lot of attention from both theoretical and ex- 
perimental sides [1-14]. This system provides a plat- 
form to study strongly correlated many-body physics in a 
highly controllable environment. The underlying interac- 
tion Hamiltonians can be engineered by diverse methods. 
This engineered system, on the one hand, can be used to 
simulate various other strongly correlated systems which 
are less controllable for achieving a better understanding 
of the involved physics, and on the other hand, can im- 
plement new model Hamiltonians which may show novel 
strongly correlated phenomena [8]. 

In this paper we describe a technique to control ultra- 
cold atoms in multi-band optical lattices for investigation 
of a class of strongly correlated physics related to the fa- 
mous Kondo problem. The Kondo problem arose from 
study of the magnetic impurities in metals [15]. Most 
of the latest investigations in this direction concentrate 
on the study of the lattice version of the Kondo model 
which is believed to be important for understanding the 
behavior of heavy-fermion and high-T c superconducting 
compounds [15,16]. Here, we show that the Kondo in- 
teraction naturally arises in the ultracold system. By 
applying spatially periodic Raman pulses, we can con- 
trol the atomic population in each band of the optical 
lattice. We show that with an integer filling of the low- 
est band, the underlying dynamics of this system is de- 
scribed by the Kondo-Hubbard lattice model, which is 
a hybridization of the Hubbard and the Kondo lattice 
models. The interaction parameters in this model can 
be well tuned by controlling the depth of the optical lat- 
tice, and this controlled realization could shed new light 
on understanding of this complicated theoretical model. 
We can realize both bosonic and fermionic versions of 
the Kondo-Hubbard lattice model, and in the case of 
bosonic atoms, we solve the model through a mean-field 
approximation, and the result shows a clear phase tran- 
sition from the Ferromagnetic superfluid to the Kondo- 



singlet insulator, arising from competition between the 
condensate-mediated magnetic interaction and the local 
Kondo interaction. 

In condensed matter systems, the effective Kondo lat- 
tice model comes from perturbation of the Anderson lat- 
tice model in some parameter region [15,16]. A recent 
work has proposed an interesting scheme to realize the 
Anderson lattice model with fermionic atoms in a de- 
signed optical superlattice [12], which could lead to an ef- 
fective Kondo model under precise control of some inter- 
action parameters. In our approach, the Kondo interac- 
tion comes from a completely different origin. In a multi- 
band optical lattice, surprisingly, the Kondo interaction 
can be derived directly from the atomic collision interac- 
tion, even in the case that the latter is spin- independent. 
As here the Kondo interaction is a non-perturbative ef- 
fect, it is much stronger compared with other approaches 
and should be ready for observation with the available 
experimental technology. In particular, we note that in- 
triguing experiments have been reported on control of 
dilute (non-interacting) atoms in multi-band optical lat- 
tices [17-19], and extension of that ability should allow 
for demonstration of the model proposed here. We also 
note the recent interesting proposals on implementation 
of the impurity Kondo model with ultracold atoms in 
different contexts [10,11]. 

The ultracold atoms considered here can be either 
bosonic or fermionic. These atoms are first loaded into 
the lowest band of a 3-dimensional optical lattice, which 
is formed by standing- wave laser beams [1]. The optical 
potential barrier is high enough so that the atoms in the 
lowest band experience no tunneling. We then transfer 
part of the atomic population to an upper band, and the 
remaining atomic population in the lowest band is con- 
trolled to be one per each lattice site. We will shown 
later how to achieve this through optical control. The 
atoms in the upper band undergo tunneling. If we ne- 
glect the atomic collision interaction, the non-interacting 
Hamiltonian for this system takes the following diagonal 
form 



k.a 



k,<r 



(1) 



where the bosonic or fermionic annihilation operators &ko- 
and fl^ correspond respectively to the atoms in the low- 
est or the upper bands, with the Bloch wave vector k and 
the spin component a. For simplicity, we consider the 
atoms with two relevant Zeeman sublevels, for which the 
effective spin component a takes two values \ or J.. For 
example, the atoms could be the fermionic Li 6 [F = 1/2) 
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or the bosonic Rb 87 or Na 23 (F = 1) while the Zeeman 
sublevel \Mp = 0) is made de-coupled by raising its en- 
ergy [20]. The lowest band energy e; is k- independent as 
the tunneling for this band is negligible. 

The atomic collision interaction in free space is de- 
scribed by the Hamiltonian 

H T = X s I (r) *t , (r) ^ (r) (r) ) (2) 

a,a> J 

where A s is the interaction parameter and ^ a (r) is the 
field operator with the spin a. For simplicity we have as- 
sumed that the collision interaction is spin-independent 
as the spin-dependent collision terms are typically 
smaller by orders of magnitudes [21]. In a two-band op- 
tical lattice considered here, we should expand the field 
operator as * CT (r) = £\ [b ia wi (r - rj) + a ia w u (r - rj)], 
where wi (r — r<) (w u (r — r,)) are the Wannier functions 
centered on the site i for the lowest (upper) band, and 
Pi* = E^P^e-^/VN {(3 = a,b) with N being the 
number of lattice sites. Substituting this expansion into 
Hj, we can derive the expression for the total Hamilto- 
nian H = Hf + Hi . Under a number of approximations 
specified below, the Hamiltonian H in the rotating frame 
has the form 

H = ^ eu^a^a^a + Uh'^rii (m - 1) 

k,<x i 

(3) 

i 

where v = for bosons and v = 1 for fcrmions. In 
Eq. (3), the number operator rij = X) CT a L a io-' tne 
three components of the spin operator {(3 = a, b) 

are defined by sf = (f}\^ - ^/3 a ) /2, of = 

(4r#a + ^At) /2, and o\ = -i - /^& T ) A 

and the coefficients = A s J |iu u (r — r^) | 4 <i 3 r, u c = 
4A S J \wi (r — ri)\ 2 \w u (r — r^) | 2 d 3 r (assuming that the 
Wannier functions are normalized). The shifted energy 
e«k = e«k — e« 5 where e u is the average energy in the 
upper band. If the band energy e„k takes the form 
of e u k = — 2i (cos k x a c + cos k y a c + cos k z a c ) (a c is the 
lattice constant for the cubic lattice) under the tight- 
binding approximation, the first term of the Hamiltonian 
(3) has the more familiar form — j) <j{ a \<j a j<?+H.c.), 
where t characterizes the nearest-neighbor tunneling rate 
and stands for all the neighboring sites. In de- 

riving Eq. (3), we have also made the following ap- 
proximations or assumptions: (i) the band gap A = 
e u — ei is assumed to be much larger than the coeffi- 
cients t,u^,u c which justifies the rotating-wave approx- 
imation; (ii) we keep only the on-site collision interac- 
tion terms for both the lowest and the upper bands; (iii) 
we have assumed one atom per each site for the low- 
est band; (iv) the Hamiltonian has been transferred to 
the rotating frame by dropping the free-energy terms 



3 blh„ + [e u + (2 + (-If) Uc /4] £ lCT a\ G a la . All 
these approximations or assumptions will be justified 
from our later numerical calculation of the lattice struc- 
ture. 

The Hamiltonian (3) represents a hybridization of 
the Kondo lattice model and the Hubbard model. For 
bosonic atoms, the Kondo interaction (the last term of 
Eq. (3)) is antiferromagnetic with the repulsive inter- 
action (A s > 0), and ferromagnetic with the attractive 
interaction (A s < 0), while the reverse is the case for 
fermionic atoms due to the different commutation rela- 
tion. Typically, the Kondo interaction is comparable in 
magnitudes with the Hubbard interaction. However, the 
ratio Uh /u c can be tuned by controlling the depth of the 
optical lattice and/or by going to different upper bands. 
With change of the optical lattice depth, the tunneling 
energy t can also be adjusted sensitively (see Fig. ID). 
So, a variety of interesting physics could be associated 
with the Hamiltonian (3) as we change the interaction 
types or tune the ratios of different kinds of interaction 
strengths. Such an extent of the controllability, if experi- 
mentally realized, would lead to significant breakthrough 
in investigation of the Kondo physics. 

Before studying the ground-state properties of the 
Kondo-Hubbard lattice model (3), first we would like to 
show how to control the atomic population in each band 
as we have assumed, and to justify the approximations 
made in our derivation of the Hamiltonian (3). For this 
purpose, we need to know the exact band structure of 
the optical lattice. Though the lowest band can be cal- 
culated easily through the harmonic approximation to 
the potential wells, this approximation is in general not 
valid for upper bands, and we find it more convenient to 
numerically solve the exact band structure. 

The trapping potential from standing wave lasers has 
the form V (r) /Er = Vq (sin 2 kox + sin 2 kay + sin 2 k z), 
where ko is the laser wave vector and we take the atomic 
recoil energy Er = Ti'k^jlm (m is the mass of the 
atom) as the energy unit so that the potential bar- 
rier Vb is dimensionless. The band structure is deter- 
mined by solving the one-particle Schrodinger equation 
-V 2 $/fc 2 + [V (r) /E R ] $ = E$, where E is the band 
energy in the unit of Er. This equation can be re- 
duced to the 1-dimensional Schrodinger equation by sep- 
arating variables $ (r) = $ x (a;) $ y (y) & z (z), and we 
can expand the wave function in each direction as su- 
perposition of the plane waves, for instance, <fr x (x) — 
J2 n a n,fc x e^ 2nfco+fcx ^, where k x is the Bloch wave vec- 
tor in x direction. We take the superposition coefficients 
fln,fc x as the variational parameters to numerically mini- 
mize the band energy E = E x + E y + E z . In this way, we 
find the exact band structure. With a typical potential 
barrier Vo = 30, the band energy E x is shown in Fig. 
1A for the lowest five bands. From the solution of the 
variational parameters a n ,k x , we can also construct the 
Bloch wave function Uk x (x) — ^2 n a n ,k x e 1 ' 2nkoX and the 
Wannier function w x (x — Xj) = J2 k Uk x e lkx ^ x ~ Xi>> / '\/N 
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in x direction (the same expression holds in other direc- 
tions). The Wannier function w x (x) is shown in Fig. IE 
for the lowest three bands. 




FIG. 1. (A) The band structure of the optical lattice with 
the potential barrier Vb = 30. The energy E x (in the unit 
of Er) of the lowest five bands is shown as a function of 
the normalized Bloch wave vector k x /2ko. (B) The enlarged 
structure of the 3rd band (the solid curve) compared with 
the fit from the cosine function (the dashed curve). (C) The 
illustration of the energy-selective Raman pulses which trans- 
fer one atom to the third band within the same internal state 
\g). (D) The widths AE X (proportional to the tunneling rate) 
of the lowest five bands shown as functions of the potential 
barrier Vo. (E) The magnitudes of the Wannier functions for 
the lowest (the solid curve), the second (the dashed curve), 
and the third (the dotted curve) bands. 

To manipulate the atomic population in each band, 
we propose to use spatially periodic Raman pulses to 
transfer atoms from the lowest band to a specified up- 
per band as shown in Fig. 1C. Note that a spatially 
homogeneous Raman pulse can not change the atoms' 
external state as the Wannier functions for different 
bands are orthogonal to each other. We can simply 
apply standing wave Raman pulses with the same pe- 
riod as the optical lattice. For instance, two Raman 
beams of the form of cos kox (the Raman Rabi frequency 
oc cos 2 k$x) with a frequency difference matching 
the band gap A = E X 3 — E x \ will transfer the atoms 
from the lowest band to the third band in the x di- 
rection. Note that the atoms are still in the lowest 
band for the y, z directions, so they only tunnel along 
the x direction and the resulting model is essentially 1- 
dimensional although we have a 3- dimensional lattice. 
We can get higher-dimensional Kondo-Hubbard lattice 
models by exciting the atoms to upper bands in several 
directions. For instance, two Raman beams of the forms 
of cos kox, cos koy respectively ( f2/j oc cos k$x cos k^y) 
will transfer the atoms to upper bands in both x and 
y directions, and intensity superpositions of the Raman 



pulses with Rabi frequencies dm oc cos ko (x + y) cos k$z, 
dm oc cos kg (x — y) cos kgz respectively will transfer the 
atoms to upper bands in all the three dimensions (Note 
that Clm + dm oc cos kox cos k^y cos k^z). 

We make use of the collision shift of the band energy 
to control the atom number in the lowest band to be one 
per each lattice site. Assume that we start with a Mott 
insulator state with all the atoms in the lowest band, and 
the average filling number of the lattice is between 1 and 
2. In this case, some lattice sites have two atoms while 
others have one. We apply the Raman pulses to transfer 
one atom to the upper band only for the lattice sites with 
two atoms (see Fig. 1C). This is possible because the 
collision energy shift depends on which bands the atoms 
are. One can choose a right frequency difference for the 
Raman beams so that only the desired two-photon tran- 
sition is resonant. All the other transitions are detuned 
by the difference in the collision energy shifts, which is 
typically about a few kHz [1]. So the Raman beams with 
the two-photon Rabi frequency significantly smaller than 
kHz will achieve the desired distribution of the atomic 
population. 

All the approximations made in the derivation of the 
Hamiltonian (3) are well justified by our calculation of 
the lattice structure: firstly, from Fig. 1A, the band gap 
A s» 18Er if we choose the third band as the upper 
band, which is much larger than the relevant energies 
t,Uh,u c (typically around Er). Secondly, as shown in 
Fig. IE, the Wannier functions are well localized even for 
the upper bands. For the third band with Vb = 30, the 
nearest neighbor collision rate is only 0.2% of the on-site 
collision rate, which can be safely neglected. The nearest 
neighbor tunneling has the dominant contribution to the 
band width as shown in Fig. IB (e„k is approximated by 
the cosine form). Finally, the tunneling rate (the width) 
for each band can be sensitively detuned through control 
of the potential barrier Vq as shown in Fig. ID. 

The Kondo-Hubbard lattice model shows interesting 
quantum phase transition arising from competition of 
different types of interactions. Here, as an example, 
we investigate the bosonic Kondo-Hubbard lattice model 
with u c > as it is easier to be realized experimentally 
and not encountered yet in the literature to the best of 
our knowledge. We solve the model through a mean- 
field approximation by assuming that the ground-state of 
the Hamiltonian (3) is not entangled for different lattice 
sites. This approximation corresponds to a generaliza- 
tion of the popular Guzwiller ansatz for the bosonic case 
[22,3], which usually gives good results in particular in 
the strong interaction region. 

With this mean-field approximation, the ground state 
energy of the Hamiltonian (3) per each site is given by 

Ei = -AE U ^2 l( a i<r)| 2 + u h (n? - Ui) + u c (s ib ■ s ia ) , 

a 

where AE U is the half band width, which is z c \t\ (z c is 
the coordination number) if only the nearest neighbor 
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tunneling is taken into account. The energy Ei can be 
numerically minimized by assuming a variational form for 
the ground state |*. t > = E mBff l ra t)„, \ m i) ai \ a )b> 
where a =|, J,, , \m{) a ^ are number states of the 

modes a^a^, and c mn<7 are variational parameters. The 
minimization is subject to the constraints ("F^ Y&i) = 1 
and (vE'il ni = n, where n is the average filling num- 
ber of the upper band. The properties of this system can 
be determined from the variational ground state. The 
particularly interesting quantities include the superfluid 

magnitude defined as S sup = J^a \{ a ta)\ I \J a\ a a icr , and 
the magnetizations of the lowest and the upper bands 
given respectively by M b = \ (s ib )\, M a = \(si a )\. 





A A E . 1 « c B A E , / u „ 

FIG. 2. (A) The superfluid magnitude (the solid curve), 
the lowest band magnetization (the dotted curve), and the 
upper band magnetization (the dashed curve) shown as func- 
tions of the energy ratio AE u /u c with the upper-band mean 
filling number n = 1/2. For this calculation, we have taken 
Auh/uc = 1.52/(4 x 0.84) = 0.45, corresponding to the real 
value with the potential barrier Vq = 30. (B) The same as 
Fig. 2A except that the mean filling number n = 1. 

Figure 2 shows the calculation results for the quanti- 
ties defined above with the filling number n = 1/2 and 
n = 1, respectively. At the integer filling, clearly there 
is a phase transition at the point AE u /u c ~ 0.5 (see 
Fig. 2B). If the tunneling rate (characterized by AE U ) 
is below this threshold value, each site is occupied by 
two atoms at two different bands, forming a local Kondo 
signet (lTDba — UT)&a) t° minimize the energy of the 
anti- ferromagnetic coupling Sj& ■ s ja . So there are no su- 
perfluid magnitude and no magnetizations in both bands. 
As soon as the tunneling rate across this threshold, both 
the superfluid magnitude and the magnetizations quickly 
go up. In this case, the atoms in the upper band actually 
form a ferromagnetic condensate to minimize the kinetic 
(tunneling) energy, while the atoms in the lowest band 
are magnetized in the reverse direction to minimize the 
energy sa ■ s ja . There are long-range correlations both 
in magnetization and in the superfluid phase. However, 
at the non-integer filling, all the quantities change con- 
tinuously (See Fig. 2B). As soon as the tunneling rate 
is nonzero, the holes in the upper band can move in the 
lattice and there is no energy penalty for that. To max- 
imize the tunneling effect (the kinetic energy), it is bet- 
ter that all the unpaired atoms in the lowest band are 
polarized along the same direction. So the lowest band 
magnetization and the superfluid magnitude starts from 



a significant non-zero value, while the upper band mag- 
netization needs to gradually increase form zero as most 
of the atoms in this band are still paired in the Kondo 
signets. 

At the end of the paper, we would like to briefly discuss 
how to observe the above phenomena in experiments. 
The superfluid magnitude can be measured through ob- 
serving the atomic interference as in the experiment [1]. 
The atoms in different bands can be separated through 
the Landau-Zener tunneling by accelerating the lattice 
with an appropriate speed [17,18,13]. To detect magne- 
tization in each band, the atomic spin states can be mea- 
sured through a Stern-Gerlach experiment or through 
spin-dependent light absorption. Note that the energy 
scale in our model is characterized by t, Uh, u c , which are 
typically about a few kHz. The temperature achieved al- 
ready in the experiment [1] should allow demonstration 
of this model. 

This work was supported by the Michigan start-up 
fund and by the FOCUS center. 
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